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Abstract 

The nuclear physics relevant to the electric dipole moment (EDM) of the deuteron is addressed. 
The general operator structure of the P- and T-odd nucleon-nucleon interaction is discussed and 
applied to the two-body contributions of the deuteron EDM, which can be calculated in terms of 
P- and T-odd meson-nucleon coupling constants with only small model dependence. The one-body 
contributions, the EDMs of the proton and the neutron, are evaluated within the same framework. 
Although the total theoretical uncertainties are sizable, we conclude that, compared to the neutron, 
the deuteron EDM is competitive in terms of sensitivity to CP violation, and complementary with 
respect to the microscopic sources of CP violation that can be probed. 
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I. INTRODUCTION 



In the field of particle physics an atomic physics quantity plays a privileged role: the 
electric dipole moment (EDM), which violates parity (P) conservation and time reversal 
(T, or equivalently CP) invariance. The Standard Model predicts values for EDMs that are 
much too small to be detected in the foreseeable future, and hence a nonzero EDM is an 
unambiguous signal of a new source of CP violation [l|, Ql- 

Over the years, many experiments have searched with increasing precision for a nonzero 
EDM. The most sensitive experiments measure the precession frequency of the spin for 
neutral systems, such as the neutron or an atom, in a strong electric field. The limit on the 

n 

EDM of the neutron, in particular, has been improved spectacularly over the years p|. The 
most precise value obtained so far is dn = (—1.0 ± 3.6) x 10^^^ e-cm j^. New experiments 
using high-density ultracold neutron sources are being set up to target the precision level of 
10-27 to 10-28 e-cm at LANL (LANSCE), PSL ILL, and Munich (FRM-II). 

Limits on the EDMs of charged particles p|, such as the electron and the proton, has 
so far been derived from experiments with selected neutral atoms (and molecules). The 
best limit for an EDM has been obtained for the ^^^Hg mercury atom j^, for which dug = 
(—1.06 ± 0.49 ± 0.40) X 10-28 e-cm was measured. In such a closed-shell atom with paired 
electron spins, the EDM of the atom arises mainly from the EDM of unpaired nucleons 
and from T-odd interactions within the nucleus. For this type of experiments with neutral 
atoms, the EDM signal is severely suppressed due to the screening of the applied external 
electric field bv the atomic electrons, a general result known in the literature as Schiff's 
theorem E, JQ- 

Recently, a new highly sensitive method has been proposed to directly measure the EDMs 
of charged particles, such as the muon or ions, in a magnetic storage ring Q, [ll|- The 
method evades the suppression of the EDM signal due to Schiff's theorem, and works for 
systems with a small magnetic anomaly. An experiment using this method has been proposed 
to measure the EDM of the deuteron at the 10-^^ e-cm level 3|- From a theoretical point 
of view, the deuteron is especially attractive, because it is the simplest system in which the 
P-odd, T-odd (f'T) nucleon-nucleon (AW) interaction contributes to the EDM. Moreover, 
the deuteron properties are well understood {l^], so reliable and precise calculations are 
possible. 
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It is our goal in this paper to address the nuclear physics part of the deuteron EDM 
calculation, and to compare the result to the EDM of the neutron (and proton) evaluated 
within the same framework. The framework is designed so that our results for the nucleon 
EDM and the f't NN interaction would be suitable, when combined with a realistic strong 
NN interaction, as a starting point for a microscopic calculation of the EDM of more complex 
systems, such as the mercury atom. 

This paper is organized as follows. In Section II we construct the general operator 
structure of the P- and T-odd NN interaction and from it derive the potential in terms of 
strong and f't meson-nucleon coupling constants. In Section III, we used this f't potential 
in combination with modern NN potential models to evaluate the two-body (polarization 
and exchange) contributions to the deuteron EDM. The one-body contributions, i.e. the 
EDMs of the proton and the neutron, are calculated within the same framework. Finally, the 
results are discussed and conclusions are drawn in Section IV. In the Appendix we discuss 
and evaluate the, also P- and T-odd, magnetic quadrupole moment (MQM) of the deuteron. 



II. P- AND T-ODD TWO-NUCLEON INTERACTION 

By contracting two Dirac bilinear covariants containing at most one derivative, the P- 
odd, T-odd, and C-even (hence still CPT-even) contact AW interaction can be constructed 
from {i) the scalar-pseudoscalar (S-PS) combination, AW x Ni'j^N, and {ii) the vector- 
pseudovector (V-PV) combination, N'y^N x W9^75W The tensor-pseudotensor (T- 
PT) combination, Wa'^'^Wx Wa^^ 75 W, also qualifies these symmetry considerations, however, 
it is equivalent to the S-PS one by a Fierz transformation. 

Using the nonrelativistic (NR) reduction and writing out the isospin structure explic- 
itly, the most general form of the low-energy, P- and T-odd {^T), zero-range (ZR) WW 
interaction, H^^^\ can be expressed, in configuration space, as 

H^^f^ = I (ci + di) cr_ + (c2 + (^2) Ti ■ T2 cr_ + (cs + ds) a + (04 + ^4) rf cr+ 

+ (C5 + rfs) (3 rl - Ti ■ T2) I ■ V (r) , (1) 

where cr± = <Ti ± <T2 and t± = Ti ± T2.^ Terms involving the isospin operator i (ti x T2)^, 

^ We note that this most general NR form containing five independent isospin-spin operators has already 
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even though they conserve charge, are ruled out since they are C-odd. The dimensionful 
couphng constants Cj and di {i = 1,...,5) each correspond to a unique isospin-spin-spatial 
operator in the S-PS and V-PV parts, respectively. These constants are the quantities 
that experiments such as nuclear EDM measurements can hopefully constrain, and thus 
predictions from different models of CP violation could be tested against with. 

At first sight, it seems that the introduction of the (ij's is redundant because the V-PV 
form has exactly the same NR limit as its S-PS counterpart, a point which has been made in 
Ref. Therefore, as long as one works strictly in the context of contact interactions, e.g. 
"pionless" effective field theory, only five coupling constants are needed to fit to experiments. 
However, there are several reasons to justify this larger set, especially when one goes beyond 
the ZR limit with several energy scales involved. 

First, when one tries to connect the experimental constraints to underlying theoretical 
models, it is still necessary to make the distinction between the S-PS and V-PV sectors. 
Because of different nucleon dynamics involved, the separation and comparison of these two 
sectors are of interest. 

Second, if one wants to keep the pions, as the lightest mesons, explicitl y an d model the 
long-range (LR) interaction through one-pion exchange (see e.g. Refs. Q,Q,Q|), a scale 
separation defined by the pion mass naturally occurs. In this case, one has in total eight 
independent coupling constants: five in the ZR potential which is a result of integrating out 
all degrees of freedom except the pions, and three /'T^ pion-nucleon coupling constants (see 
below, Eq. (jH)) which describe the LR potential.^ This possible scale difference between the 
S-PS and V-PV sectors is not manifest in the ZR limit. 

The third and more practical reason is that we are going to adopt a "hybrid" approach 
for the AW dynamics which takes advantage of existing high-quality strong NN potentials 
and use perturbation theory based on operators constructed in the spirit of effective field 
theory (EFT). In such a framework, it is necessary to smear out the contact interactions. 
The physical guideline is to take the delta function as a limit of the mass^-weighted Yukawa 
function, m^3^j,.(r) = ml e~"^''^ / (Anr) , when the exchanged boson is taken to be extremely 

been pointed out in Ref. Q|. 
^ The three ttAW couplings were first pointed out by Barton However, since the concern then was 
parity violation, these couplings were only picked up later when interest in nuclear CP violation built up. 
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massive: 



lim ml = lim F.T. 



ml 



5^3) (r), (2) 

where "F.T." stands for Fourier transform. As suggested above, allowing different mass 
scales for the S-PS and V-PV sectors then leads to the most general Hf^ in terms of ten 
independent operators. 

Although the choices of the mass parameters for the Yukawa functions are arbitrary in the 
sense of fitting the coupling constants, the mass spectrum of low-lying mesons provides an 
intuitive choice and suggests a connection between H^f thus constructed and the one-meson 
exchange scheme. Besides the one-pion exchange (J^ = 0", m^, = 140 MeV) often adopted 
in the literature, the contribution from r] (J^ = 0~, = 550 MeV) and from p and 
uj (J^ = 1~, mp^i^ = 770,780 MeV) |3| have also been considered in various works. We 
will show that a one-meson exchange scheme containing vr, r/, p, and oo produces the same 
general operator structure as the ZR scheme. (The isoscalar-scalar meson e or "cr", with a 
^'t coupling of type Ni'j^aN, leads to the same operator structure as the rj meson, and its 
contribution would be effectively subsumed in the coupling G^\) 

The strong and /'T^ meson-nucleon interaction Lagrangian densities, Cs and Cf,f, are^ 

'Cs = gTTNNNi-f^T ■ 7zN 

+gr,NNNi-i^7]N 

-9.NNN{l^^ - i^a^'q^y^N , (3) 



and 



+N{-g(^)r^ + -gWr^r^)N 



+ £'V^^^)a^^g,75iV . (4) 



JV 



^ The choice of pseudoscalar coupHng for the pion field in Eq. ^ is traditional in the EDM literature. In 
order to have manifest chiral symmetry, pseudovector (derivative) coupling is of course preferred. However, 
the results for the two-body contributions and for the leading one-body contribution (the chiral logarithm) 
would be equivalent. 
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The gxNN^s are the strong XNN coupling constants for which we will adopt the values: 
g^r^N = 13.07 HQ, 9vr^N = 2.24 Q, 9pnn = 2.75 [2^ and q^^^n = 8.25.^ The (^ij^'s are 
the ones with the superscript i = 0, 1, 2 denoting the corresponding isospin content.^ 
Xv and Xs are the ratios of the tensor to vector coupling constant for p and u respectively; 
when vector-meson dominance (VMD) is assumed, they are equal to the electromagnetic 
counterparts, i.e. Ky = 3.70 and Ks = —0.12. The tensor structure N a^" q^'-^^N in Eq. (jij), 
where = Pu —pl, is equivalent to the PV structure iV9^75A^ by a Gordon decomposition. 

Evaluating all one-meson exchange diagrams with one strong and one vertex, the NR 
potential, Hff, is found to be 



+rf cr , ■ V 



+ (3r,V,^ - n ■ T,) <T . V(G(^) y^{r) - Of 3^,(r))} 



(5) 



where is defined as the product of a strong coupling constant Qkmn and its associated 
ft one for instance, G^^^ = Qt^mnQ^^ ■ 

One sees that the general operator structure in Eq. (|T}, based only on symmetry consid- 
erations, is fully reproduced by the one-meson exchange scheme containing the lowest-lying 
pseudoscalar and vector mesons in both isovector (vr and p) and isoscalar {rj and u) sectors. 
The ten coupling constants in Eq. find their counterparts in the ten ft meson-nucleon 
coupling constants. Eq. ^ has the advantage that it not only has the most general operator 
structure, but it also provides a link to the meson-exchange picture which provides some 
insight. We finally note that one-kaon exchange does not contribute to the strangeness- 
conserving NN interaction. 



We use the prediction [2^ gt nn = ^9pnn to infer guiNN from gpt^^/'in = 0.6 given in Ref. 

We use the Bjorken-Drell metric and special attention should be paid to the definition of 75 = ^ ^ ^ ^ 

and any coupHng constant associated with it. Relative to the PauH metric a sign difference due to 75 
should he.^kept in mind. 

In Ref. the fit r]NN coupling only contains an isoscalar part, so it does not contribute to the isovector 
Hfj-. However, this isovector piece, which gives a different linear combination from the pion contribution, 
is needed in order to render the <t_ and rl (t_|_ operators independent. 
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III. DEUTERON EDM 



Because the f'f interaction induces a small P-wave admixture to the deuteron wave 
function, it leads to a nonvanishing matrix element of the charge dipole operator. In addition, 
since the proton and the neutron also have an EDM, a disentanglement of one- and two-body 
contributions, 

djy = dg) + d!§ , (6) 

is necessary to make contact to the underlying f't physics. In the following, we shall use 
the f't NN interaction Hf>^ constructed in the previous section to calculate dg'*. We will 
use the same meson-exchange picture as a guideline to give an estimate of . The final 
result for dx> can then be expressed in terms of the meson-nucleon coupling constants. 
EDMs are expressed in units of e-fm for the remainder of the paper. 



A. Two-Body Contributions 

For the two-body part, the dominant contribution comes from the polarization effect: In 
leading order in the perturbation, it is the matrix element of the charge dipole operator 
evaluated between the unperturbed deuteron state \V) (mainly ^Si-wave with some 6% 
^Di-wave) and the admixed P-wave component viz. 

dt''' = \j\mr'-er\\V), (7) 

where r = ri — r2 and "||" denotes the reduced matrix element. Because the charge dipole 
operator conserves the total spin, \V) has to be the ^Pi state. The isospin and spin selection 
rules then dictate that only the operator rf; cr+ in Hff can induce such an admixture to the 
deuteron. 

In order to examine the model dependence of the matrix element, the numerical calcu- 
lation is performed with three high-quality local potential models: Argonne Wig (Afig) j^l, 
and the Nijmegen models Reid93 and Nijm II (2^. The results 

= 1.43 X 10-2 G« + 1.59 X 10-3 aj') + 6.25 x 10-"^ G« - 5.96 x 10-^ , (8a) 
= 1.45 X 10-2 G^^^ + 1.68 X 10-3 G^^^ + 6.83 x 10-^ - 6.53 x 10"^ , (8b) 
= 1.47 X 10-2 Gt^) + 1.72 X 10-3 gJ^) + 7.50 x 10-^ GJ,^) - 7.19 x 10"^ G[^) , (8c) 
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for Avis, Reid93, and Nijm II, respectively, show a relatively model-independent pattern. 
Judging from the coefficients for the different mesons, pion exchange dominates the result. 
The much smaller sensitivity of d^^ ~ d^"''^ to heavy-meson exchanges guarantees that pion- 
exchange is a good approximation here (this may not be true for d^\ a point we address 
below). 

The slight difference in the results for these models can be attributed to their softness 
at the intermediate range where the deuteron wave function (which agrees well for these 
potential models) has most of the overlap with the Yukawa functions. Fig. [T] compares the 
effective potential Ve//(r) = Vs{r) + L{L + l)/(m^rr^) of these models in the ^Pi channel 
(L = 1), which determines the radial behavior of the ^'f admixture in the inhomogeneous 
Schrodinger equation 

{T + V}P^)\V) = H,,\V). (9) 

Among these three models, Nijm II is the softest one within the range of about 0.3 ~ 
1.0 fm, so it gives the largest result, while At>i8, the hardest one, gives the smallest result. 
As the heavy-meson exchange is very sensitive to the wave function at short range, its 
model dependence is more apparent compared to the pion-exchange case. Our result for 
the coefficient of is consistent with two earlier predictions: 0.010-0.026 obtained by 
Avishai 3|, who used strong potential models from before the 70s, and 0.019 obtained by 
Khriplovich and Korkin who assumed the zero-range approximation for the deuteron 
and a free ^Pi wave function. Their number can be considered as an upper bound. 

The meson-exchange effects, in the form of two-body exchange charges, give contributions 
of the form 

dg") = (^{V\\ J d^xp^^\x)x\\V) + {V\\ J d^xp^^\x)x\\V)^ , (10) 

where the first term corresponds to adding the normal (P- and T-even) exchange charge 
p*^^-* to Eq. (0), and the 1^7^ exchange p*^^\ induced by H^f^ is included via the second term. 
Compared with the one-body charge, which is 0(1), p^^-* can be ignored since it gives a 
correction of (9(l/m^r)^ (see e.g. Ref. |2i]). On the other hand, since p*-^^ can be as large 
as (9(l/m^), and its contribution is evaluated within unperturbed deuteron wave functions, 
its significance should be investigated. 

As indicated by the dominance of pion exchange observed above, and also in view of the 
suppression of heavy-meson exchange currents found in the study of the (P-odd, T-even) 
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Figure 1: The comparison of three different effective strong potentials in the ^Pi channel. 

deuteron anapole moment j^^l, the consideration of the pion sector is sufficient for the two- 
body exchange effects. Attaching a photon to every possible line in the one-pion exchange 
diagram vv^hich leads to H^J , the exchange charge can then to C(l/m^) be identified, in 
configuration space, as 



((1 + Ks) [Gf T,-T, + G« + (3 ~ n ■ r,)) 



+(1 + K,) (G(°) + + 2 r^)) a, ■ V. 5^'\x - n) y^{r) 
+ (1^2), (11) 



-(tt) / _ \ 



Ami 



(12) 



v\fhere the pair term refers to the diagram in which the photon couples to an intermediate 
nucleon-antinucleon pair and the mesonic term refers to the diagram in which the photon 
couples to the meson in flight; r = jri— r2|, rxi(2) = \x — ri{2)\- Numerically, the contribution 
of these diagrams to the deuteron EDM is found to be 



— 



9.40 X 10"^ GW - 5.28 x 10-^G1°) . 



(13) 



Compared with d^^°^\ this constitutes only a few-percent correction. 



Combining the results for and we obtain for the two-body contribution to the 



deuteron EDM, in terms of f't couplings, 
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(14) 



N 



N 



(xy - -0 

p,co % 



N 



N 



Figure 2: The two-body contribution to dx) arising from the first-order f^'f p- and u-TT'y couphngs. 

with an error estimated as less than 5%. 

Besides the usual exchange effects in which one of the meson-nucleon couplings is P- 
and T-odd, another class of diagrams involving a |^7^ photon coupling to the exchanged 
mesons can also contribute. Since pseudoscalar mesons cannot have such a coupling 
to photons, the candidates in our current framework are pn'y, un'y, and pp7 vertices. 
Assuming these |^7^ couplings are of the same order of magnitude, one can expect a smaller 
contribution from the pp7 vertex, because the p meson is much more massive and has a 
smaller strong coupling to nucleons than the pion. Therefore, in order to estimate the size 
of this type of contributions we evaluate the diagrams based on J^T^ pn^y and un'y vertices 
shown in Fig. [2l 

Expressing the pny and ujtt'j Lagrangian densities as 



C 



(pTT-y) 



2m 



(15) 
(16) 



where two new /'T^ coupling constants (^^^^ and Qi^.,^^ are introduced, the associated exchange 
charges are, in configuration space. 



Pmesonici^'i ^1) ^2) 



Pmesonici'^'i ^1' ^2) 



n ■ T2 (V ■ V2)(<T2 ■ V2) ypir,i)y^ir,2) + (1 ^ 2^7) 

;v ■ V2)(tT2 ■ V2) yu.{r.i)y.{r^2) + (1^2). (is) 



The numerical calculation, using the Afig potential, gives an EDM contribution of about 
2.3 X 10^'^ {dpn-y — duj-K-y)- Slucc the coefficient is two orders of magnitude smaller than the 
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Figure 3: Hadronic loop diagrams which contribute to the nucleon EDM. 



leading coefficient oi g^^^ in Eq. (fTll) . we shall ignore these mesonic f'f effects for the rest 
of this work. 



B. One-Body Contributions 

The total one-body contribution to the deuteron EDM is simply the sum of the proton 
and neutron EDMs, i.e. 

=dp + dn. (19) 

Our goal in this section is to evaluate dp and dn in a manner consistent with the framework 
used for the f'^t NN interaction. 

The nucleon EDM has a wide variety of sources such as the QCD 9 term, quark EDMs 
and chromo-EDMs (CEDMs), Weinberg three-gluon operator, and four-quark contact in- 
teractions, therefore, its evaluation requires good knowledge of nonperturbative dynamics 
of confined quarks, which is still not available. A commonly used method of estimate 
is to evaluate the hadronic loop diagrams, in which meson and baryon degrees of free- 
dom are used to describe the dynamics and the dependence on the |^7^ mechanisms at the 
quark-gluon level is subsumed in the J^T^ meson-nucleon coupling constants. This approach 
has been applied extensively to the neutron EDM in various contexts (see, for example, 

Refs. 0, Q, E, Q H Q, Q, Q, Q)- 

Here we apply it to both the proton and the 
neutron EDM, with the inclusion of vector mesons. 

The loop diagrams containing a virtual pseudoscalar meson are classified as in Fig. El^a) 
and Fig. in^b), where the photon couples to the charged pseudoscalar meson in the former 
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and to the intermediate nucleon in the latter case. Defining the hadronic loop contribution 
to the nucleon EDM as 

(20) 



Ahad) 



the results for the corresponding diagrams are^ 



1 t"" 



0. 



(21) 
(22) 



<5f = (3Gf + GW)jW-(3/^,Gf + «:,G«)jW 

+ (G(°) + G«)j;'') - («:,GW + «:.G«)T?) , (23) 
5^ = (-Gf + G« + 4Gi2))jf) - (-«:vGW + + 4«:,Gf )J(-) 

+(Gf + G(i))X;'') - («:,Gf + KsG^^W . (24) 

The three distinct loop integrals involving an i-type pseudoscalar meson, 

correspond to the cases where the photon couples to the pseudoscalar meson, the nucleon 

Dirac, and the nucleon Pauli form factor, respectively. They are evaluated as 



t 



= -1 - (1 - x^) Inx, + (3 - x,2) ^(x^) 
> -Inxi - 1 + ^ Xi + xf In + O(x^) , 



-^1 



1_1 2, 



^1 (1 - 2^^^) -^(^ 



1 



1 



^(S) 



Y^ + g^i - Y^a;^ (l + 3xj Inxi-Xi (1 + -^^ 
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x^<^i 16 4 
1 



V4s + s2 



x^ Inxj + 0(x: 
2-s 



tan ^ 



+ tan 



(25) 

(26) 

(27) 
(28) 



_y/4s + s\ 

where Xj = rrii/m^. Since both meson and nucleon form factors were taken to be constant 
off the mass shell, and since form factors fall off as the square of the four-momentum transfer 
increases, these results should be viewed as an upper bound 

From Eqs. (pK|l -|?7jl. one observes that only Xq'^^ has a non-analytic term, i.e. InXjr, in 
the chiral limit, — > 0. The mathematical reason is that Fig. Efa) contains more pion 



'' Kaon loops can also contribute 



411 |. and be easily added to our results. 
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propagators than Fig. Of b), which is responsible for the infrared divergence in the soft-pion 
limit 2J|- Therefore, the contribution to the nucleon EDM involving chiral logarithms is 
purely isovector, 

This implies that the deuteron EDM receives no one-body contribution from loop diagrams 
involving tt and rj mesons in the chiral limit. Furthermore, when the neutron EDM is 
considered, the constant terms in To and Ji exactly cancel, as has been pointed out in 
Ref. 0|. However, this is not true for the proton. 

Because chiral symmetry is explicitly broken by the pion mass, it is interesting to compare 
the chiral logarithm with other, analytic, terms when realistic parameters are used. For 
example, taking Xj^ = 140/940 in Eq. (f^ . one gets 2q^^ ~ 1.19, which is about 40% smaller 
than — ln(m^/miv) — 1.90. This sizable difference typically sets the scale for the theoretical 
uncertainty. The same conclusion can be drawn from the work by Barton and White 
who, motivated by the success of sideways dispersion relations for the nucleon Pauli form 



factors 



42|, applied the same technique with the same parameters to the neutron EDM 



problem. This analysis, involving mainly the threshold pion-photoproduction amplitude, 
is actually similar to the evaluation of type (a) loop diagram with soft pious, and in fact 
produces the same chiral logarithm. Compared with the leading term which gave (i„ ~ 
0.8 X 10^^^ e-fm (a value, in our notation, glr^ —gl?^ = 1-2 x 10^^° was used as input), their 
full analysis predicted dn — 0.5 x 10^^^ e-fm, which is also some 40% smaller. While this 
may be just an accident, it does signal a potentially large theoretical uncertainty for the 
nucleon EDM. 

In order to estimate the relevance of the vector-meson degrees of freedom to the nucleon 
EDM, we consider the diagrams illustrated in Fig.[ni[c). These contributions can be roughly 
estimated by the assumption of VMD, which leads to a dispersion-theory analysis of the p° 
and u poles in the time-like region. The deuteron is only sensitive to the isoscalar sector, 
for which, in the case of the nucleon Pauli form factor, the naive VMD model works rather 
well. The vector-meson contributions to the isovector nucleon EDM should, however, also be 
added as a correction to the leading result from the pion-loop calculation, which is equivalent 
to including the 27r continuum in the dispersion-theory analysis. 

The required vector-meson-photon conversion mechanism is introduced by the Lagrangian 
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density 

C — P/^^ f'^p) -u P^"^ pH f^nl 

where the F^'^'s denote the field tensors for the photon and the p° and uj mesons; the 
constants /, =^00 and /„ = 17.05 are determined from the decay widths r,_...- = 
6.85, 0.60 MeV |42| by T^^e+e- = m^r/(3/J). Then the vector-meson contributions to 
the nucleon EDM are evaluated as 

= J^LgW + ^G™, (31) 
= -^(GT + + 7^G« . (32) 

Jp9pNN Ju)9uJNN 

Keeping in mind the caveat of a possibly large theoretical error, we nevertheless take 
a more adventurous point of view and include also the analytic terms in Xq and Xi, i.e. 
X^"^ ~ 1.19, Xf^ ~ 0.41, and x{''^ ~ 0.28, but we neglect the part from the nucleon Pauli 
form factor, so we set X^^^ = X^''^ = 0. Collecting the results from Eqs. (l2TT)-(l2ll. 
(|3T| - (j32ll . the total one-body contribution to the deuteron EDM is evaluated as 

dg^ = 2.18 X 10"=^ (3 + G-W) + 1.49 x 10"=^ (Gf + G^) 

+1.53 X 10-2 g(^) + 1.49 X 10"^ G(,°) . (33) 

In terms of the ^'t meson-nucleon coupling constants, the result is 

4^ = 0.03 ^« + 0.09 # + 0.04 ^« + 0.01 + 0{£'''^) . (34) 



IV. DISCUSSION 

Combining Eqs. (|T4l) and (|34| . we arrive at our final estimate for the deuteron EDM: 

dv = (0.20 + 0.03) ^« + 0.09 # + 0.04 ^« + 0.01 £^ . (35) 

while our results at the same time imply the following predictions for the proton and neutron 
EDMs: 

d, = -0.08(# - #) + 0.03(# + -gl'^ + 2#) + 3 x 10-^(4°) + ^W) (36) 

-f0.02(^f + + 2^f ) + 6 X 10-^(£°) + , 
d^ = 0.14(# - ) - 0.02(^f - ^« + 2^f ) + 6 X 10-^(£°) - -g2^) . (37) 
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The leading contribution to dx>, 0.20 glr^^, due to the 1^7^ AW interaction, including the 
exchange charges and calculated by using state-of-the-art wave functions, is 25% smaller than 
the result assuming the zero-range approximation j^^l which was adopted for an analysis 
on CP violation models in Ref. (4J|. The remaining contributions come from the proton 
and neutron EDMs. These terms have a sizable theoretical uncertainty, which could be as 
large as 40%. The non- vanishing dependence on gn \ which arises from including analytic 
terms in the hadronic loop calculations, sets the stage for the QCD 9 term to play a role 
in the deuteron. Using the prediction by Crewther et al. Q| that g'^^ ~ 0.027^, one gets a 
(ix) ~ 2 X 10^^(9 dependence on 9, which is about three times larger than the QCD sum-rule 

nn 

calculation |4J, |43J. The dependence on the vector-meson couplings, though suppressed 
at the two-body level, enter the final result through the nucleon EDM where it could be 
sizable. An important issue in this respect is the size of the p, u (and rj for that 
matter) coupling constants compared to those of the pion. An argument by Gudkov et 
al. suggested that these vector-meson coupling constants are less significant while a 
recent work by Pospelov based on QCD sum rules gave the "best" values for gp^^ of the same 
order of magnitude as g.„ ji^; and these two works surprisingly have opposite predictions 
about the relative importance of g-ni/g-^o. Furthermore, work on the P-odd, T-even AW 
interaction implied vector couplings at least equally important and preferably larger than 
their pseudoscalar counterparts (see e.g. Refs. {43, 13, Therefore, until consensus 

is reached, these vector-meson contributions should still be kept for maintaining a greater 
generality. 

In order to connect expression Eq. for the deuteron EDM to the underlying CP viola- 
tion, the ft meson-nucleon coupling constants have to be expressed in terms of parameters 
at particle-physics level, such as the QCD 9 term, quark EDMs and CEDMs, etc. These 
quantities have a plethora of predictions from extensions of the Standard Model. Because 
all the EDM measurements to date only resulted in upper bounds, it is a popular practice to 
use these experimental limits to derive the corresponding bounds for one particular source 
of CP violation while turning other possibilities off in an ad hoc fashion. Even though this 
simplification is legitimate to some extent, one might obtain an overconstraint by excluding 
possible cancellations between various CP-violation sources. 

The deuteron and neutron results illustrate how limits on their EDMs could be used to 
provide tight constraints on a specific model of CP violation, such as the one in Ref. Q- For 
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supersymmetric models in which the Pecci-Quinn symmetry is evoked to remove the QCD 9 
term, the quark CEDMs are the dominant contributors to the f't meson-nucleon coupling 
constants, compared to the three-gluon and four-quark operators. Therefore, all the q's 
can be expressed in terms of the (ij^'s. Using the "best" values recommended in Ref. |4fil |: 
# ~ 20 di, ~ Ad\, gf^ ~ 13.3 rf^, g^p^ ~ -8.6 rfl, - -8.6 rf^, gl^^ ~ -13.3 di, 
where = (i^ic?^,^ the deuteron and neutron EDMs can be completely expressed in terms 
of the CEDMs of the up and down quarks, viz. 

dj, = -A.Q7d''^ + 5.22dl, (38) 
dn = -OMd^a + OAgdl. (39) 

Thus, these two EDM measurements probe different linear combinations of c?^ and c?^ in this 
case. Moreover, the deuteron could be significantly more sensitive than the neutron. This 
example is clearly oversimplified, however, judging from the general expressions Eqs. ()35|) 
and (f37|l . one expects that, barring unnatural and accidental cancellations, the deuteron is 
competitive to the neutron in sensitivity to CP violation. Furthermore, the deuteron EDM 
involves different coupling constants, and hence in general will be complementary with 
respect to the information about CP violation that can be probed with the neutron. 

In conclusion, it should be realized that the theoretical uncertainties, especially in the 
results for dp and dn and hence in the one-body contribution to dx>, are significant. The 
calculation of an atomic or nuclear EDM involves a broad range of physics, including the 
problematic strong interaction at the nuclear and subnuclear scale. In this respect, it is rele- 
vant that efforts have been renewed recently to attack the neutron EDM in lattice QCD 0|. 
In general, improved treatments of the hadronic physics, which can bridge the phenomenol- 
ogy of the neutron EDM and /"T^ nuclear forces with the underlying particle physics, are of 
central interest. 
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Appendix A: DEUTERON MQM 

Besides the EDM, the NN interaction can also induce P- and T-odd electromagnetic 
moments of higher multipolarity, that is, C3, C5, and M2, M4, etc. For a spin-1 object such 
as the deuteron, a nonzero M2 magnetic quadrupole moment (MQM) is therefore another 
signature of CP violation. Approximating the nuclear electromagnetic current as purely 
one-body, i.e. ignoring the meson-exchange currents, the MQM operator can be expressed, 
in a Cartesian basis, as 
g 

M^n = {/^ [3 2 (T ■ r 6mn] + '2 q {rjn Ln + rn L^)} , (Al) 



where /i, g, and L denote the nuc 
angular momentum, respectively 



eon magnetic moment, charge (in uints of e), and orbital 
The deuteron MQM, defined by 



2 

Mr, = 2{V,J, = l\Y, M,,(z)|P, J, = 1) , (A2) 

i=l 

can then be evaluated once the ^Pi and ^Pi parity admixtures have been calculated. As- 
suming one-pion exchange only, and using the Af ig strong potential gives the numerical 
result 

Mj, = 0.051 jjs # + (0.031 fiy + 0.003) g^J^ , (A3) 

in units of e-fm^. The model-dependency is at the 1% level, similar to the EDM calculation. 
Although the isoscalar spin current leads to a rather large matrix element (in the zero- 

n 

range approximation of Ref. |aO|, it is three times the isovector one), the isoscalar magnetic 
moment renders the resulting g^^ coefficient, 0.04, smaller than the gl^'^ coefficient, 0.15, 
which is dominated by the isovector spin current from the large isovector magnetic moment. 
The orbital motion adds only a small correction to the glr'^ term through the deuteron 
D-waye component. 

While a sensitive experiment to measure appears as least as formidable as for rfx), 
it might be contemplated with deuterium atoms, because the MQM, unlike the EDM, is 
not screened by the electron. An interesting theoretical point is that, since the nucleon 
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itself has no quadrupole moment, the deuteron MQM is a rather clean probe of the f't NN 
interaction, and in particular of 
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